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A review of developments in the mathematics and methods for principal value (PV) 
integrals is presented. These topics include single-pole formulas for simple and higher- 
order PVs, simple and higher-order poles in double integrals, and products of simple poles 
in general multiple integrals. Two generalizations of the famous Poincaré-Bertrand (PB) 
theorem are studied. We then review the following topics: dispersion relations for the 
advanced, retarded, and causal Green’s functions; Titchmarsh’s theorem; applications 
of the PB theorem to two- and three-particle loop integrals; and the R and T matrix 
formalism. Also, various applications of the PV methods to nuclear physics, transport 
theory, and condensed matter physics are studied. In the appendices several methods 
for evaluating PV integrals, including the Haftel-Tabakin procedure for calculating the 
R and T matrices, are reviewed. 


1. Introduction 


In recent years there has been a resurgence of interest in applications of PV methods 
to intermediate-energy nuclear physics (see Refs. 12, 13, 16, 46, 52 and 60). Most 
of these studies focus on numerical calculations using various response functions or 
“loop integrals”, which arise from nucleon-delta-pion interactions. Such a response 
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function, involving a product of two single-particle Green’s functions, then appears 
as a self-energy correction which is inserted into an external propagator, as shown 
in Fig. 1. 


G,(o + @) 


Gpl’) 


Fig. 1. Feynman diagram for a two-propagator loop integral. The loop is imagined to be inserted, 
essentially as a self-energy correction, into an external propagator (which is illustrated in the 
diagram by the dotted line).!?»52 


If one examines the dispersion relations for a general finite-temperature loop 
integral, one finds that its real and imaginary parts are related by an expression 
that involves two PV terms in a double integral.!?5? Then, in order to simplify this 
expression, it is usually necessary to reverse the order of integrations in the double 
integral, an operation that is most elegantly handled using the PB theorem (see 
Refs. 2, 12, 15, 48, 52, 57 and 60). These manipulations are discussed in Sec. 3 of 
this paper, in which we also treat the dispersion relations for the finite-temperature 
three-particle loop integral.!? (See Fig. 2.) Then, some results of numerical cal- 
culations of two-particle loop integrals in intermediate-energy nuclear physics are 
reviewed in Sec. 5 of the paper. 

There have also been many applications of PV methods in nuclear physics,?9*1 
condensed matter physics? and transport theory.3® Moreover, the PB theorem 
was used in a study, in momentum space, of the quantum mechanical reflection 
and transmission at a potential step. Also, recently®® certain PV logarithmic 
integrals, which arise in the plasma energy of an electron gas, can be evaluated 
using a generalized dilogarithm. Such research, of course, is only a small sampling 
of all of the many applications of PV methods in physics, applied mathematics, and 
engineering. 

Thus, there has been new interest in the mathematics of PV integrals (see Refs. 2, 
11, 12, 15, 44, 47, 48, 57 and 60), and in analytic and numerical methods used to 
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Fig. 2. Feynman diagram for a three-propagator loop integral. Again, the dashed line represents 
an external propagator into which the loop is inserted. 


calculate such integrals (see Refs. 10, 11, 13, 15, 17, 29, 47, 54, 55 and 60). 
The purpose of this paper is threefold: 


(i) to review the mathematics and methods of PV integrals (Sec. 2 and Appendices 
A, B and C); 

(ii) to discuss the application of PV methods to general techniques commonly used 
in physics, mathematical physics and applied mathematics, such as dispersion 
relations and various Green’s functions (Secs. 3 and 4; Appendix B); 

(iii) to review more detailed applications of the PV formalism in nuclear physics, 
transport theory and condensed matter physics (Secs. 5-7, respectively). 


Our review of PV methods and applications is by no means exhaustive. The 
literature on this subject is simply too vast to be covered in a single review paper. 
However, we do believe that the papers referenced and the topics and examples 
presented here are representative of the most important kinds of PV studies occur- 
ring in mathematical physics. Moreover, the examples that we address here largely 
involve numerical treatment of PV integrals as opposed to earlier studies in which 
PVs were evaluated analytically, often by complex variable calculus. 

Typically, most physicists learn to manipulate and evaluate simple PV integrals 
early in their professional careers. Of particular importance are infinite-limit inte- 
grals containing a single pole, which can be evaluated by complex variable theory if 
the integrand contains fairly simple analytic functions.114” However, such integrals 
can usually be evaluated almost mechanically, without any real understanding of 
the true meaning of the PV; i.e. one wishes to know the precise interpretation of 
the equation 


Ila) =? f copie (1.1) 
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where the integration limits can either be finite or infinite, g(x) is “well-behaved” 
over the region of interest, and P denotes the PV operation. 

We believe that there are a number of misconceptions regarding the fundamental 
meaning of PV. Also, for many applications, it is crucial to understand the PV 
operation. In particular, for developing accurate numerical methods to evaluate 
PV integrals, an intimate familiarity with basic definitions is absolutely essential 
(see Refs. 10-13, 17, 29, 47, 54, 55 and 60). In Sec. 2 of this paper, we attempt to 
clarify certain problems associated with misconceptions, misinformation, or plain 
mistakes about the PV operation. 

The first common misconception about the PV is that, in an equation like (1.1), 
one is somehow evaluating a contribution from the pole. This is nonsense, of 
course, since in most cases the pole contribution is infinite. (Whether the inte- 
gral in Eq. (1.1) without the PV symbol actually diverges depends on the limits of 
integration and the precise form of g(x).) In Eq. (1.1) what one is evaluating is 
the integral over the entire range of integration, except for an infinitesimally small 
region near the pole. In fact, a fundamental definition of the PV that is strongly 
suggested in the literature is that the PV is the “convergent part” of an otherwise 
divergent integral (see Refs. 11, 15, 17, 24 and 37). While this definition may be 
almost trivial for simple, single pole formulas (like Eq. (1.1)) its extension to more 
complicated cases requires considerable justification. 

A second misconception about the mathematics of PV integrals is that the basic 
definition is limited to an isolated simple pole, as in Eq. (1.1). In fact, the PV idea 
can be usefully extended to higher-order poles!?-15; e.g. 


1) (qa) = pf a. ae DT de, n>1, (1.2) 


or to two or more PVs occurring in a multiple integral; e.g., 


1o =? [5 | gew. (1.8) 


where again the limits of integration can either be finite or infinite. Equation (1.3) 
is the type of integral, previously mentioned, to which the PB theorem and its 
generalizations apply (see Refs. 2, 12, 15, 48, 52, 57 and 60). 

One important “spinoff” of the PB theorem is a generalization of the well-known 
equation** 


— 7i6(a — a) (1.4) 


(z -a+ ic)! — 
e~0 T—-—@ 


to a double product like 


(a —at ic) (x — b + teg)7?. 
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Then, as €1,€2 — 0 we find that 
(x —a+ie,) (x — b+ te)? 
P ; P ; 
# | — ni6(x — o) K — nid(x — b) ; (1.5) 


The correct relation is given by!?*? 


(z-a +icı) (x —b+ iez)! 


P P 
e120 l za ee a) Z — mid(x — 3) + 1°6(x — a)5(x — b). 
(1.6) 


The difference between Eqs. (1.5) and (1.6) arises from the PB theorem. Equa- 
tion (1.6) was recently used in Ref. 35. 

There is one final refinement of the PV definition that needs to be mentioned. 
For finite-limit PV integrals we can take advantage of certain “quadratures” in order 
to evaluate the integrals.!> In particular, assume the existence of a quadrature for 
the indefinite integral 


(n) = g(x) 
Fi™ (x) G-a dz, (1.7) 
which in principle always exists. Then, it can be shown that 
i z 
PJ m = FM (c) — F(®){(b), (1.8) 


a result that is valid whether the point a lies on the interval between b and c, or 
outside that interval. In other words, the PV of the integral is simply the differences 
of the quadratures, obtained by ignoring the singularity and “blindly integrating 
through it”. Once again, we see that the PV is the convergent part of the integral. 

With these preliminaries, we next begin Sec. 2, containing the formal PV defini- 
tion and its generalizations. The mathematics of PV integrals is covered in Sec. 2, 
and a number of important applications are reviewed in Secs. 3-7. In Sec. 8 we give 
a brief summary of this paper. 

In Appendix C we present a remarkable property of PV integrals, virtually 
unknown, that greatly extends the number of one-dimensional PV integrals that 
can be evaluated exactly. 


2. Principal Value Integrals 

2.1. Single-pole formulas for principal value integrals 

2.1.1. General formulas 

First, we explain our terminology. A single pole expression is a single integral whose 


integrand is given by 


Ga(2) = ee n>, (2.1) 
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where z is the integration variable, a is a constant, and n is an integer. Then, a 
simple pole occurs for 
n=1, (2.2a) 


while a higher-order pole occurs for 
n>1, (2.2b) 


and we also say that the pole is of order n. We will assume that the function g(x) 
is “well-behaved” in the integral. More precisely, it is assumed that the function 
and its first n — 1 derivatives are finite and continuous over the range of integration. 
(For most cases of physical interest, g(x) and all of its derivatives will be finite and 
continuous over the range.) 

For a simple pole, the fundamental definition of the PV is!!:44 


1 er (x-a) 
(saa) =e. e oe 
while for higher-order poles, we have 
wci fet a) 
penae a ie alea) aa 


Note that Eq. (2.4) reduces to Eq. (2.3) for n = 1. For a simple pole, one also has 
the well-known relations** 


1 1 : 
pee =P (—.) F ind(x — a), (2.5) 


with the limit £ — 0 on the left-hand side understood and with (x — a), the delta 
function, given by 


E 


A Ree (x =a)? +E? 


(2.6) 


We emphasize again that each of the expressions in Eqs. (2.3)-(2.6) occurs in an 
integrand involving a single integration variable, and the range of integration can 
be finite or infinite. 

From Eq. (2.3) one can show that?1-44 


ca ee 


ne ie +f leat r 


and if a lies outside of the range of integration [b, c], then the PV integral in Eq. (2.7) 
becomes an ordinary integral. The limits b and c can either be finite or infinite. 
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Formula (2.7) shows that P can be written either inside or outside the integral sign. 
For n > 1, Eq. (2.4) in an integrand can be expressed as!* 


P f eo ate= LEE o-o] 


1 € dz gir-1) 
"Ge pi? f G-a) | Z) , (2.8) 


where 
G-=1) d(n-1-3) 3 
(a) d r-a 
$; (u )= | dzíi— sani 9) fim L 1—3) | (x — a)? al i 
dG-}) qir-1-3) 1 
= - |S g(x x) dx("-1-3) (l (2.9) 


where the last line follows if u # a. Equation (2.8) is derived by a series of integra- 
tions by parts!® and we obtain from Eq. (2.7) 


© oft) 4 Oo Suae a) 
Pir ae nA [HO - a0] 


tui) mis a +f I< (x —a) “oa are DIEE (2-10) 


Thus we have reduced the higher-order PV integral to a simple PV expression. 
Invariably, this is the technique used to evaluate all higher-order PV integrals.*112-15 
Also, in many cases the boundary terms in Eq. (2.10) vanish. We also note that 


Ane ee (2.11) 
dx*) (x ~ a) “(x — a)kHl? ` 
so that Eq. (2.9) becomes 
(yy) [Eaa du 
; (u) = maa G-I £) (2.12) 


t=u 


Thus in Eqs. (2.7) and (2.10) we have the equations for simple and higher-order 
pole PV integrals, respectively. In fact, we can consider Eq. (2.7) as a special case 
of Eq. (2.10) since, for n = 1, the boundary terms in Eq. (2.10) do not occur. There 
are two special cases of interest which we shall now consider. 


2.1.2. Infinite-limit integrals and the connection with complex variable theory 


We let b — —oo and c — +00 in Eq. (2.10). We also assume that the boundary 
terms all vanish, i.e. 


®;(t00) = 0, (2.13) 
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and that g(x) — g(z)(with z = x + iy) is an analytic function in the upper-half 
plane. Then, by the standard techniques of complex variable theory!! 4”, Eq. (2.10) 
can be expressed as 


mi d?-) 
Py... ea z C dr ~(n—1)! —1)! l dz("-1) goo) 9 ’ (2.14a) 
where a is assumed to be real. If g(z) is analytic in the lower-half plane, we have 
—ri di” 1) 
RLA G-a (ni | a) sot (2.14b) 


Equations (2.14) are the usual expressions one gets for a pole, with the PV integral 
differing from the usual residue expression! 47 by the fraction, 1/2. For a simple 
pole, we have the important expression 


hee 2p l PNE iy (2.15) 


=œ (x z a) 


which is the Hilbert transform of g(x). 

However, we should mention that there exist infinite-limit integrals, which one 
cannot easily evaluate using complex variable theory. An excellent example of such 
an integral is the Hilbert transform of a Gaussian, namely” 


1 co go? 


and e~*’ is well behaved along the entire real axis, but has nontrivial analytic 


behavior in other parts of the complex plane. However, by a standard manipulation 
of this integral,? it can be shown that 


h(a) = -e J En (2.17) 


the right-hand side of which can be expressed in terms of the error function of 
imaginary argument. We emphasize that Eq. (2.17) can be derived without any 
assumptions on the analytic behavior of e7, 

Consequently, we see that the complex variable theory definitions of the PV1147 
are not the most general or even necessarily the most useful. For many cases of 
„interest, one can easily evaluate infinite-limit, PV integrals using complex variable 
techniques. Such methods are particularly useful if the function g(z) can be simply 
separated into a function which is analytic in the upper-half plane and another 
which is analytic in the lower-half plane. When such a separation is not possible, 
other methods may have to be adopted to evaluate the integral.? (For example, see 
Ref. 39.) 
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Also, we now see again the motivation for the (n — 1)th derivative in Eq. (2.4). 
Basically, the derivative in the PV is converted to a derivative on the well-behaved 
function g(x) by a series of integrations by parts, giving us the general expression, 
Eq. (2.10), which involves only a simple pole PV integral. However, the PV of a 
higher-order pole is automatically obtained from complex variable theory,!147 in 
which the natural definition of a residue is the (n — 1)th derivative of g(z) evaluated 
at the pole. If the pole is along the real axis, we get the PV, with a factor of 1/2 
times the “usual” residue,!1*47 as in Eqs. (2.14). (The usual residue is for a pole 
entirely within the closed contour, not on its boundary.) 


2.1.3. Finite-limit integrals and quadrature relations 


In Eq. (2.8) we now specialize to the case in which b and c are finite. Consider the 
indefinite integral of the left-hand side of Eq. (2.8), namely 


_ [ _9(2) 
F(z) = f G= Phas (2.18) 


The function F(z) is then a “quadrature” of the integral. If g(x) is sufficiently 
simple, F,(z) can usually be expressed in “closed form”. However, from the funda- 
mental theorem of integral calculus,4 F,(z) always exists in principle. 

If we perform a series of integrations by parts on the right-hand side of Eq. (2.18), 
we obtain?” 


(-)" n—1 ; di-)) din—1-3)( 1) 
F,(x) = CET ry eos) oe 


1 dz dir) 
‘mes go) : a 
Equations (2.9), (2.10) and (2.19) suggest the following relation 
° glz) 
Pf de = F,(c) - Fa (b), 2.20 
Pinte = Ful) - Fuld) (2.20) 


so that PV is simply the difference of the quadratures. In Ref. 15, Eq. (2.20) is 
rigorously derived for both simple and higher-order poles. 

Also, Eq. (2.20) shows that the PV is obtained by ignoring the singularity, i.e. 
“blindly integrating through the pole”. As we will show more formally in the next 
subsection, this is because the PV is really the convergent part of the integral. 
Thus the next subsection is basically the rigorous justification of the derivation of 
Eq. (2.20). 

There are a number of simple PV integrals that can be readily evaluated by 
quadratures, e.g.1° 


° dr c 
Pf 7 = In|; |, (2.21a) 
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€ dæ ted 
a z= + re (2.21b) 
dz 1 1 
, oe OT OR Oe) 


Also, it is sometimes convenient to shift the pole in Eq. (2.20) to the origin,!® 
i.e. 


F, (c) — F, (b) = P f OY gola) ay (2.22) 
b—a T 
where 
Ja(x) = o(x +). (2.23) 


2.1.4. The PV as the convergent part of a divergent integral 


Consider a general integral of the form 


(a) 
Ia(b,c) = > i dz ES (2.24) 


where n refers to the highest- or pole in the integrand. (Also, b and/or c can be 
finite or infinite.) Note that g (x) may contain additional powers of (x — a), e.g. 


| ķi ean i cond p L giy [e7 (x — a) + e7], 


so that g(x) = = e~°* (x — a) and gS? (2) = zeme, 
Now in Eq. (2.24) we can expand each g?(z) in a Taylor series, namely 


BOET AO ENEE | 


= (E)E — a! E 


k=0 


dg} (x) 
dx L i 


dé 
EO s)| . (2.25) 


=a 


Then, a sufficient condition for the convergence of the integral Ia (b, c) is the follow- 
ing 


n dk j 
2 EX ES =0; k=0,1,...,n— 1. (2.26) 
j= 


However, note that each separate term in the expansion on the right-hand side of 
Eq. (2.24) is divergent. 
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Moreover, from Eq. (2.11) we see that 


— (yt fe a) (1 Vw 
I,(b,c) = Saf dt ze (£). (2.27) 


If Eqs. (2.26) are satisfied, we can integrate by parts n — 1 times the right-hand side 
of Eq. (2.27) to obtain 


cor 1 N n-1 (a) qir-1-k) 1 @r=C 
Ta(b,¢) J! L (Sa 93 (2) X Fey ma] 
j=l k=1 x=b 
N 
© ded) (o 
J =) ena) deers e: (2.28) 
I= 


Since [,(b,c) converges, there is no loss of generality in taking the PV of the 
right-hand side of Eq. (2.28), so that 


)* N n-1 
I, (b,c) a DS —)* 
j=1 k=1 
dik-1) (2) q(n—1-k) 1 2=e 
TEDI (2) e Dn 
dz d(-) (a) 
(n i= wand $? [ (x — a) aj dx) 95 (x )- (2.29) 
Then, from Eqs. (2.9) and (2.10) we find 
Ane Pf (a 2.30) 
( ) 2 E (z-a) _ ayr 99 ( ye ( 


We thus obtain the important result that a convergent integral, consisting of terms 
each of which is separately divergent, can be evaluated by summing the PVs of 
each of its terms.!!!5 Thus, the PV can be interpreted as the convergent part of 
the integral. 

We also assume that if the limits |b| and |c| are infinite, the appropriate boundary 
terms in Eq. (2.29) will vanish. There are many examples of both finite-limit and 
infinite-limit integrals that can be readily evaluated using Eq. (2.30) (e.g. Refs. 11 
and 15) 


R(t) Jêl 71 
al9,Cc oa wT el el 2.31 
I,(b,c) = [aw [els 2 Béla (2.31) 
where 6 and a are constants and 
R(x) = [2? + ye + 8]. (2.32) 
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The entire integral converges since Eq. (2.26) is satisfied, but each of the three 
terms on the right-hand side of Eq. (2.31) separately diverges. Also the range of 
integration is assumed to exclude the zeros of R(x). It can then be shown that! 


R 2 
I(b,c) = _R(2) + In|2R(x) + 22 + y| — In 2/6|R(x) + 26° + y2 
x 2|6| x 
jól y g=c 
ra zl as (2.33) 


For infinite limit integrals (b — —0oo,c — +00), an additional rearrangement of 
the integrand is desirable, namely separations into terms that are convergent in the 
upper- and lower-half planes. Then Eqs. (2.14) can be used to evaluate the various 
terms making up the integral, e.g. 


i sin z — z COS £ d=3 f sin z — xcosx d 
A z3 ae a Sm r3 j 

B 1° et? — eit (et? + e`?) 

=. \( in )-2 z? Je 

1 og isli 1 1 Pa isf i 1 
-Peker Lel(e-2) 


(2.34) 


Each term in integrand of the first line of Eq. (2.34) has a pole of order 2, but 
it is clear that this singularity disappears when the two terms are combined. Note 
too that, after we make the separation into terms proportional to e+** and e~** 
(which are analytic in the upper- and lower-half planes, respectively), there appear 
terms having poles of order 3. However, each of the final four terms may be easily 
evaluated using Eqs. (2.14), giving the final result. : 

We also remark that the material presented in this section is almost intuitively 
obvious to most physicists and applied mathematicians. In fact, most scientists 
instinctively use the approach outlined here when evaluating a convergent integral 
consisting of a sum of separately divergent integrals. However, a formal derivation of 
this method seems to be missing from the standard literature, and the justification 
is particularly important for cases in which one encounters higher-order poles. 

Also, related discussions can be found in the mathematics literature concern- 
ing the rigorous theory of generalized functions?*3” and of Hadamard finite-part 
integrals.1” However, the language is different, and the application in generalized 
function theory (which is usually expressed in terms of “regularization” of divergent 
integrals) may not be entirely equivalent to our approach. 

In Appendix A we summarize all of the many methods (see Refs. 17, 13, 15, 29, 
54, 55 and 60) that have been developed for numerically evaluating PV integrals. 
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In all such numerical methods, we make use of the property that the PV is the con- 
vergence part of the integral. By a subtraction method, we explicitly eliminate the 
singularity in the integral, so that the convergent PV part can be readily evaluated. 


2.2. Formulas for multiple integrals; the Poincaré-Bertrand theorem 
and its generalizations 


2.2.1. Simple poles in double integrals 


There are many examples where two PVs occur in a double integral (see Refs. 12, 15, 
35, 52 and 60). In manipulating such integrals, one must be careful about reversing 
the order of integration. In particular, the PB theorem? !?,4*-5” states the following: 


O TEAGI 


P z 
=f af ap ogien- ie (23) 


where f(x,y) is well-behaved on the real axis. The last term on the right-hand side 
of Eq. (2.35) is the result of reversing the order of integration. 
This theorem is also valid for finite-limit integrals!548-57, e.g. 


c P ff. P 
[ezl dy ED) 


af nf up o P 
=f av f EA 


=r? f(u,u) foc — u)O(u — b) + O(u — c)O(b — u) |, (2.36) 
and O(w) is the Heaviside step function 


1 fw >0 
Oe) = e if w <0. 20) 


The limits of integration b and c are constants, independent of x and y. The final 
term involving the step functions in Eq. (2.36) assures us that z = u or y = u lie 
within the ranges of integration. If u is outside the ranges of integration, then both 
PVs become ordinary integrals, and the order of integration can be reversed without 
the need of the extra term. Other finite-limit expressions, analogous to Eq. (2.36), 
can easily be constructed (e.g., for the case in which the z and y limits are different). 
In Ref. 15 an example of a finite-limit, double PV integral is presented. In particular, 
it is shown analytically that the last term in Eq. (2.36) disappears if the range of 
integration does not include u. Other demonstrations of the theorem for finite-limit 
integrals can be performed!® in numerical evaluations. (See Appendix A of this 
paper.) The derivation’!:** showing the equivalence between Eqs. (2.3) and (2.7) 
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(thereby verifying the validity of interchanging the PV operation and the integral 
symbol) can be extended to the double integrals in Eqs. (2.35) and (2.36). 

The general rule that one must remember in reversing the order of integration 
in a double integral involving two PVs is that if one or both of the PVs disappears, 
the extra term does not occur. The elimination of a PV may occur due to the limits 
of integration, as we have seen in Eq. (2.36). However, a PV can also be eliminated 
if f(x,y) in Eqs. (2.35) and (2.36) is proportional to factors of (x — u)”, (y — u)” 
or (x — y)” (with n a positive integer > 1). In what follows we will always assume 
that both PVs are present in the double integral. 

There are two useful identities that occur in this formalism. First, we have!248 


EE ee eRe ee | oa a So 
rora raes eal: ee 


Equation (2.38) can formally be justified using the fundamental PV definition, 
Eq. (2.3) and then taking a limiting procedure. However, Eq. (2.38) is almost 
obvious. Clearly, it is valid if there are no PVs present (ie. if z Æ u,y # u, 
and y # x). Recall from Sec. 2.2 that, when one takes a PV, one is eliminating 
the divergence due to the pole. Thus, basically Eq. (2.38) is valid because of the 
“regularization” property inherent in all PV expressions.?*37 

The second important identity is a generalization of the familiar Eq. (2.5). If 
there are two poles present in an integrand, one can show that,!249 


[(x — uy — 181€1)(%@ — ug — isse2)]~* + ims,6(x — u) 


— 
€1,€2-0 | T — U2 


z +inszô(z — ua)| + r?s(x — u,)6(@ — ua), (2.39) 
— U2 


where s; = +1. The last term in Eq. (2.39) arises from the last (extra) term present 
in the PB expressions (2.35) and (2.36). It is clear that in using Eq. (2.39) in an 
integral, a PV must not be eliminated for any reason; e.g., the range of integration 
must include both u, and ue. 

Examples of the usefulness of the PB will be discussed in Secs. 4-6. However, 
there is a simple mathematical example which we now present. Let h(x) be the 
Hilbert transform of g(x). Then, from Eq. (2.15) we have 


Lf P 
ae . 2.40 
na) == f _ valu) (2.40) 

Next, construct the function 
ss he fe, OP 
aw) = 5 f de- Mo) 
Le fo oe Pe SD 

anya nae 2.41 
= fe [aw (2.41) 
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which from Eq. (2.35) becomes 


Z 1 P 
Gu) = 1 g(y)dy he dr —— a 7 Gan — g(u). (2.42) 


It can be shown that 


oF P P 
J EE TE = 0 (u real), (2.43) 


whereby Eq. (2.42) becomes 


ĝlu) = —g(u) = f° do—_W(2). (2.44) 


Thus, we obtain the well-known result!:?247 that if h is the Hilbert transform of g, 
then g is (minus) the Hilbert transform of h. This result is usually proven by com- 
plex variable theory by letting z — x + iy and requiring that Y(t? (z) = h(z) +ig(z) 
be an analytic function in the upper-half plane. Here the result has been established 
without any assumptions on analyticity. Nevertheless, we know from a fundamental 
theorem due to Titchmarsh*® that there is an intimate connection between analytic- 
ity and the existence of Hilbert transforms. We shall discuss Titchmarsh’s theorem 
in much more detail in Sec. 3. 

2.2.2. Higher-order poles in double integrals 


The PB theorem can also be generalized to higher-order poles. It has been shown 
that}? 


a ad EPA > if ae g yten) 
-faa toy te ha 
-emen ("7 Vara ten] 
(2.45) 


where 


@) E w (2.46) 


Note that Eq. (2.45) reduces to Eq. (2.35) for n = m = 1, for which the only 
contribution to the sum on the right-hand side of Eq. (2.45) is for r = 0. 
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In the derivation}? of Eq. (2.45), the following partial fraction expression was 
used 


P P P Ser 
(z — u)” (z =y)” eo a 

| P P 

@ =u y 


mt+r-1 


LPs ey e 
r'=0 


(m+r-1-r')! 
P P 
1—1 2.47 
e E E a 
However, it is clear that the left-hand side of Eq. (2.47) is invariant under the 
interchange of u and y and of n and m. Thus, we find that!4 


P P (-)” N (n+r-1)! 
ere e 


r=0 r! 
AS 
G—y Gay 
n+r—1 
n+r - „(m+n-r —2)! 
nema E 
P p 
x onra (2.48) 


However, because of the symmetry of the left-hand side of Eq. (2.47), yet a third 
(manifestly symmetric) expression can be derived,!* namely 


ae a APE. Ne pp eel) P P 
(x —u)™ (x — y)™ = (m—1)! on r! (z — u)" (y — u)" 
(-)” l (n+r-1)! P P 
i (n-1)! >D (x =y) (u—y)rtr ` (2.49) 


Note the high degree of symmetry present in Eq. (2.49) as compared to Eqs. (2.47) 
and (2.48). All three expressions are equivalent! and can be used to derive vari- 
ations of the PB theorem for higher-order poles. Also, the remarks regarding the 
derivation and validity of Eq. (2.38) also apply to Eqs. (2.47), (2.48) and (2.49). 


2.2.3. Simple poles in general multiple integrals 


The final generalization of the PB theory is for multiple PV integrals which contain 
only simple poles.!? We define the following functions: 
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P P P 
ARP (w) ~ J (zı — zef (z2 — zef (£3 — a” 


xX xX [oe Jn(£1,£2,£3,..., 2n) (2.50) 


(tn — 21) 


BY) (u) = [den [dena f des f des 
P P 


fn(21,22,°*+ En), (2.51) 


and 


* Gi — u) (En — 21) 


and for simplicity of notation we do not put the limits on the integrals since the 
formulas that we present pertain either to finite-limit or to infinite-limit integrations. 
(As in Subsec. 2.2.1, it is assumed that no PV is removed in any of the integrals.) 
The function f,(21,22,...,;2n) is assumed to be well-behaved over all of the n 
domains of integration. 

By comparing the definitions in Eqs. (2.50) and (2.51) with Eq. (2.35), it is 
clear that the natural generation of the PB theorem to multiple integrals can be 
expressed as 


AP (u) = BP (u) + DY (u), (2.52) 


where DÑ Vw) is the extra term resulting from exchanging the various orders of 
integration in going from AY Yu) to BY (u). Of course, for n = 2 we have from 
Eq. (2.35) 


DP (u) = -r° fa(u, u), (2.53) 


and, after extensive manipulations, one finds that!? 


DP) = 2? | 


(z 2 po [ote J f3(u, x, u) = fa(u, uÙ, x) . (2.54) 


Now, a general expression for DP (u) can be written but it is rather unwieldy and 
not very useful. Therefore, we shall instead give a rather practical recursion relation 
for finding DY Vu) for relatively small values of n (say, n < 5). 

Suppose that we have a function of n + 1 variables, namely 


Jn+1(21, £2, ee Eni Entia 


and we eliminate the last variable by defining another function 
P 
In(T1, Tr, wage Zn) = dtng fing (21, T23-- -3 Ln, Tn+1) : (2.55) 
(n41 = Tn) 


Now assume that DY \(u) is known for gn(£1,£2,...,£n). Then, it can be shown 
that}? 
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DY) (u) = DO) (u) +(x? | dan f den- vf des 


P 
X | =r fn. (U, £2, £3,- .., Ln, U 
Lee at +1(U, £2, £3 ) 


EEEE E E 
a (i = u) Mele = 2) 


= 


+ tay tet) . (2.56) 


We note that in the last two lines of Eq. (2.56) the first and last variable positions in 
fn+1 are taken up, first, by u and then by zı. Also, we do not integrate over either 
1 OF Ln41. This asymmetry in coordinates results from the special link between 
zı and £n in Eqs. (2.50) and (2.51). As an example, consider a function f3(z, y, z) 
for which we define 


92(2,y) = i Z dzfa(2,y,2), (2.57) 

and from Eq. (2.56) we find that 
DP (u) = DI (u) +r? f dy | 
From Eqs. (2.53) and (2.57), this becomes 
DY (u) = =r? J 


+a faiy[ punt tono] es 


u-y 


f(u,y,u) + = fon] . (2.58) 


u-y 


ya qifun Us z) 


which is Eq. (2.54). Then, since we have the general structure of DY Vu), we can 
find DY Vu) and higher-order DYP functions. Of course, in practice one is limited to 
rather small values of n (say, n < 5); otherwise, Eq. (2.56) becomes very difficult to 
use. Note also that the asymmetry in the coordinate positions in Eq. (2.54), again 
resulting from the asymmetry inherent in the original equations (2.50) and (2.51). 

We remark too that the P’s appearing in the second and third lines of Eq. (2.56) 
have the following meaning!” 


P 


AOP n 
z Å K”, 2.60 
I- (2 — 2) A (x-xi) ( ) 
where 
-1 

KP = | [ħa] . (2.61) 

j=1 

j+i 


Thus, a PV involving a product of factors can be expressed as a sum of the PVs 
of each of the individual functions. These remarks, coupled with the results in 
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Eqs. (2.47)—(2.49) imply that we can always express a complicated PV (involving 
a product containing poles of arbitrary order) in a partial fraction expression. In 
fact, a general expansion for a PV of the form 


aa eee 
Ta (2 — ti)" 


is given in Ref. 14. In this case, the expansion is in terms of individual higher-order 
PVs such as 


(N, ni > 1) 


P 
ENA i21). 
(x — t;)™ (m: 2 1) 
Equation (2.60) yields another interesting identity. It is well known that 
E E ai Gane = a real. (2.62) 
Thus, from Eqs. ou and r we find that 
=0 (x; all real) . (2.63) 
One other important Boie should be stated, namely 
© dr 
P per = 0, (2.64) 


which can be proven by a partial fraction expansion and a limiting procedure. 

Finally, it should be mentioned that in Ref. 12 other generalizations of the 
PB theorem to multiple integrals were considered (i.e. generalizations other than 
(2.50)-(2.52)). It was concluded that such generalizations are difficult to construct 
and appear to be uninteresting and useless. On the other hand, the generalization 
given here has many applications, an example of which is given in Sec. 3. 


3. Dispersion Relations 


3.1. Dispersion relations for advanced, retarded, and causal 
Green’s functions 


It is very useful to have a relation between the real and imaginary parts of a Green’s 
function. For example, in the case of a finite-temperature causal Green’s function 
(or propagator), there is an important dispersion relation?’ 


Re [G(w -if a! = Ace nw! (3.1) 


where the Green’s function, G(w), is in the energy representation (E = fw) and 


nw) = { tanh [(w — 4)/2kpT] for bosons 


coth [(w — u)/2kgT} for fermions; 2) 


pis the chemical potential, kg is the Boltzmann constant, and T is the temperature. 


852 K. T. R. Davies et al. 


In the zero temperature limit, we have 
nw) z sign (w — p). (3.3) 


Then, from Eqs. (2.40) and (2.44) we find that 


fiiC@inwy === °° he aa IGT (3.4) 


T (w — w) 
For the advanced and retarded Green’s functions, the above formulas also apply, 


except for Eq. (3.2), which becomes 


(3.5) 


ye +1 for the retarded function 
MW) = _1 for the advanced function. 


Sometimes it is useful to resort to the time representation using the Fourier 
transform 


oo 
G(t) = J dw e~™'*G(w). (3.6) 
—o0 
The inverse of the Fourier transform (3.6) is 


Gwi= = f 7 dt eG). (3.7) 


3.2. Titchmarsh’s theorem 


The theorem® states that, for f;(z) and f2(x) Lebesgue square integrable functions 
for all values of z, the following properties imply one another: 


(a) F(z) = fi(z) + ife(z),z = 2 + iy, is an analytic function in the upper-half 


plane, with 


J | F(z) |? dz finite for y > 0. (3.8) 


(b) fı(x) and f2(x) are Hilbert transform of one another, with 


n=? pee’), (3.9) 


=œ (z' zi z) 


(c) The Fourier transform of F(z) is proportional to a step function. Specifi- 
cally, identifying z with an energy or frequency representation, x > w, we 
find that in the Fourier transform (time) representation 


F(t) = f 7 dw e~*”* F(w) (3.10a) 
= Q(t)G(t). (3.10b) 
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This last property is the causal condition. The validity of any one of these three 
conditions implies the other two. Thus, we see very clearly from this theorem the 
equivalence of the existence of a Hilbert transform, analyticity, and causality. 

Also, we mention that for a function that is analytic in the lower-half plane, the 
Hilbert transform alters sign in Eq. (3.9) and Eq. (3.10b) is replaced by 


F(t) = O(-t)G(t). (3.11) 


It is instructive to give a partial proof of Titchmarsh’s theorem since some 
results of this proof will be useful in further work.4! We do not attempt to establish 
the connection between condition (a) and conditions (b) and (c), so we make no 
assumption about analyticity or square integrability. However, we will demonstrate 
that conditions (b) and (c) imply one another. 

In the proofs that follow, the reader may notice that the order of integration is 
reversed for double integrals containing a PV. This does not contradict the results 
of Sec. 2.2 since the integrals in question contain only one (not two) PV. In fact, in 
integrals containing products of Schwartz distributions? (e.g., PVs, delta functions, 
or step functions), the order of integrations can always be reversed provided that 
only one PV is present. (That is, in a double integral containing two delta functions 
or a delta function and a PV, the order of integrations can be freely exchanged.) 

First, assume Eq. (3.10b), where F(t) has the Fourier transform 


FEA i i iwt 
Flu) = > f _ Pitta. (3.12) 
Next construct the function 
o0 1 
a=2f mee Telas (3.13) 
S Ww) 
which becomes 
eit 
Hojs zm | I G(t)dt P i n aa" . (3.14) 


However, from complex variable theory we know that 
et S 
P — yy! — = = iwt : 
w E T mi [O(t) — O(—t)] et, (3.15) 
which gives 
1 ita ; 
I P ; iwt 
(w) = Im l; f G(t)e at 


=n l 5 re iwt 
= Im if Fi(tje t| 


= Re [F(w)}], (3.16) 


the desired dispersion relation. Thus, condition (c) implies condition (b). 
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Now assume the existence of the dispersion relation 


Re[F(w)] = È 2S aw as EEA, (3.17) 


-w) ’ 


so that 
F(w) = Re[F(w)] + i Im [F(w)] 
1 1 P ! 1 
z. / dis [5 tine -.)] Im[F(w’)], (318) 


which from Eq. (2.5) becomes 
_1 f° ,, Im[F’)| 
F(w) = = [iw aol: (3.19) 


Now, the Fourier transform of F(w) is given by 


ev iwt 
zf du’ Im [F(w" fz = (3.20) 
(w —w ie) 
Then, since we have 
o0 —iwt ae 
lim —* dw = +2rie™™tO(t), (3.21) 


e0t J œ W — w — ie 


Eq. (3.20) becomes 
F(t) = 2i0(t) T dw'e™™ t Im F(w’), (3.22) 


which is proportional to the required O(t). 

Thus, we see that the existence of a Hilbert transform implies causality, and 
vice versa, with no particular assumptions about the analytic properties or square 
integrability of F(w). The complete proof of the theorem®® establishes the final 
connections, linking analyticity, causality, and the existence of Hilbert transforms. 

A useful identity also follows from Eq. (3.22) of Ref. 41. First, we note that 


F(w) = Re[F(w)] + i Im [F(w)] 
= fiw) +ife(w), (3.23) 
from which we obtain 
F(t) = falt) + ife(t). (3.24) 


Of course, f,(w) and fo(w) are real functions by definition, but f(t) and f2(t) 
are clearly complex. Then, from Eqs. (3.10a) and (3.22)-(3.24), we find that*! 


filt) = ~ife(t) + 120(t) f(t) 
= ifa(t) [-O(t) - O(—t) + 20(8)] (3.25) 
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or 
ft) = i [O(t) - O©(—t)] fo(t) (3.26) 
or 
folt) = ~i [O(t) - O(—t)] A(t), (3.27) 
where we have used the identities 
1 = O(t) + O(-t) (3.28) 
and 
@(t)@(—t) =0. (3.29) 


Also, from Eqs. (3.24), (3.26), (3.27) and (3.28) we obtain*! 


F(t) = O(t) fo(t) (3.30a) 
= 20(t) fı (t). (3.30b) 


Thus, we see that, in either the w or t representations, there is only one inde- 
pendent function (either fı or f2). This behavior, of course, arises because of the 
existence of the Hilbert transform, Eq. (3.9). Also, as before we make no assump- 
tions about analyticity in deriving Eqs. (3.26), (3.27), and (3.30). 


3.3. Application of the Poincaré—Bertrand theorem to response functions 
or two-particle loop integrals 


The response function,” or loop integral,!?52:° is proportional to the energy- 
representation integral 


Teisi / T duo! Ga (w +02")Gelw!). (3.31) 


(See Fig. 1 for the appropriate Feynman diagram.) In Eq. (3.31) Ga and Gg are 
the propagators for particles a and (, respectively. In the time representation, 
Eq. (3.6), Eq. (3.31) becomes 


Maa(t) = iGa(t)Ga(—t) (3.32a) 


or 
Maa(t - t’) = iGalt — t')Ga(t' -t), (3.32b) 


which accounts for the loop occurring in the response function. 

We assume now that œ and @ pertain to Green’s functions that satisfy the 
dispersion relations (3.1) and (3.4). For causal functions, no(w) and ng(w) are 
given by Eq. (3.2), while for retarded and advanced propagators we use Eq. (3.5). 


E E E ERE 
Then, the imaginary part of Eq. (3.31) is given by 
maslo) = f” dw {Re[Ga(w + o) RelGo(w") 
— Im[Ga(w + w')JIm[Ga(w')] }. (3.33) 
From Eq. (3.1), we obtain 
a du! Re[Ga(w + w')] Re [Gs(w’)] 
=a au! [7 diva g — y In [Gala] 


x Mala) [dog — eae (3.34) 


Now, in Eq. (3.34) we can interchange the w’ and wa integrations since only one 
PV is present,” and from Eq. (2.35) we find that 


i dw’ Re [Galw + w’)] Re [Gg(w’)] 
z 5 T dwa Im [Ga(wa)] nala) f dwgIm [Ga(wa)] ng(we) 


meee P P 
x [as (wa —w'—w) (wg — w’) 
J dwa Im [Ga(wa)] Na(Wa)Im [Gg(we — w)| ng(wa —w). (3.35) 


From Eq. (2.63), the first term in Eq. (3.35) vanishes and after redefining variables, 
Eq. (3.33) becomes 


Im [Meg(w)] = T dwa S dwg 5(Wa — w — wg) 
x [na(wa)ne(we) — 1] Im [Ga(wa)] Im[Ga(wg)]. (3-36) 


From Eqs. (3.1) and (3.31), we obtain 


Re [Ia(w)] = -> A 7 die i ae 


P 


X — 
(Wa — w — we) 


x [Ne(Wa) — na(we)) (3.37) 


Im [Ga (wa)] Im [Golwg )] 
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and it is clear from Kgs. (3.36) and (3.37) that the following dispersion relation is 
valid 


Re[Mlaa(w)] == f- di! ln Ha") al") (3.38) 


where 


nelwg) ~ Malwa) (3.39) 


Tap lwa ~ wg) = Na (Wo )Na (we) =1 


If a and ĝ pertain to finite-temperature causal Green’s functions, then from 
Eqs. (3.2) and (3.39) we obtain! 


tanh [(w — Ha + ue /2Kk BT] 
(ie if a and 8 are particles with like statistics , 
eB) = | coth [(w — pa + a)/2kBT] 
if a and @ are particles with unlike statistics. 


(3.40) 


This result states that if both particles are bosons or both are fermions, then the 
loop integral behaves as a boson. On the other hand, if one particle is a boson and 
the other is a fermion, the loop behaves as a fermion. In the zero-temperature limit, 
we have 


Nap(w) z sign (w — Ha + up). (3.41) 
If the propagators are retarded and advanced Green’s functions, we find that 


0, if both propagators are 

retarded or advanced functions, 
+1, for a retarded, 8 advanced, 
—1, for a advanced, ĝ retarded. 


Nop(w) = (3.42) 


Equation (3.42) is a consequence of the physics which demands that, in the 
zero-temperature limit, only particle-hole combinations are allowed in the response 
function. In fact, this result can be seen more generally from Eqs. (3.36) and (3.37) 
if we take!? T — 0. In that limit, particle-particle and hole-hole contributions 
identically vanish. Of course, for finite temperatures, the particle-hole boundaries 
are smeared, and one cannot make an unambiguous identification of a particle-hole 
pair. 


3.4. The Poincaré-—Bertrand theorem applied to three-particle 
loop integrals 


It is of interest to consider the generalization of the PB theorem to triple PV in- 
tegrals, for which one uses Eqs. (2.50)-(2.52), for n = 3, and Eq. (2.54). Such 
analyses then give much insight into how to generalize the results to multiple in- 
tegrals of arbitrary order. Two special cases of interest were treated in Ref. 12: 
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(i) convolutions with respect to time and (ii) triple-propagator loop integrals. In 
this section we will discuss only the latter. 
In Fig. 2, we display the triple-propagator loop integral, which is given by 


Lapy(w) = T dwa I dwgGa(wa)Gg(we)G,(we +wp +w). (3.43) 


Then, from Eqs. (3.6) and (3.7), it can be shown that in the time representation we 
have 


Lapy(t) = Ga(—t)Ga(—)G,(t), (3.44) 


and the direction of the arrows in Fig. 2 clearly correspond to the time arguments 
in Eq. (3.44). 
From Eqs. (2.50)-(2.52), (2.54) and (3.1), we find that}? 


P co o0 co 
Re [Las (w)] = =f dwa | dws f dwy 


Im [Ga(wa)] Tm [Ga(wg)] Im [G(y) 
(w + We +wg-— wy) 
X [Na(wa)na (wa), (wy) — Nnalwa) — ngolwg) + ny(wy)]. (3.45) 


Also, one can show that 


Im[Lagy(w)] = o dwa Í. dwg 
x Im[Ga(wa)] Im[Gg(wg)] Im[G (wa + wg + w)] 
x [ne (we)ny (wa + wg +w) +n (wa + wg +w)nalwa) 
-= No(we)na(w_e) — 1). (3.46) 


From Eqs. (3.45) and (3.46), we obtain the three-particle loop integral dispersion 
relation!” 


p i Im [Lap (w')] 
Re [Lapy(w)] = = J È dus" Napy(w') on (3.47) 
where 
Nosy (Wa, wp, wy) 
aby (Wy — Wa — Wg) = aAA, 3.48 
Nesp (wy a) Daou lata.) ( ) 
with 


Napy (wa WB, Wy ) = Na (wa )ne (wg ny (wy) — Na (wa ) — ng (wa) + Hy (wy) (3.49) 


and 


Dopy(We, We, wy) =1+% (wa )na(we) = ne (w_e)n,y (wy) — Na(We Im (wy) . (3.50) 
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Then, for a finite-temperature causal function it can be shown from Eqs. (3.2) and 
(3.48)—(3.50) that 


tanh{w — (Hy — Ha — He)/2kBT} 
ye for three bosons or two fermions and one boson, 
“lee coth{w — (Hy — Ha — Ha)/2kBT} 
for three fermions or one fermion and two bosons. 


(3.51) 


By comparing Eqs. (3.2), (3.40) and (3.51) we see a very satisfying physical 
result. The statistics for the loop depend upon the “statistical combination” of 
the particles making up the loop. For example, two fermions combine to form a 
boson, while the composite of three fermions is a fermion. Then, in the dispersion 
relations Eqs. (3.38) and (3.47) the appropriate 7 function is a tanh (with a suitable 
argument) for boson-like loops and a coth for fermion-like loops. Clearly, these 
results can be easily generalized to more complicated types of loop integrals. 


4. The R and T Matrix Formalism in Momentum Space 


We consider scattering solutions of the time-independent Schrödinger equation 
(T+V)V¥p = EŬp, (4.1) 


and T and V are the kinetic and potential energy operators, respectively. The 
solutions of Eq. (4.1) can be expressed in the Lippman-Schwinger formë? 


Ye =O24+GeV¥r, (4.2) 
where ®z is a solution of the equation 
(T — E)$z = 0, (4.3) 


and Gz is a Green’s function involving the inverse of the operator (T — E). In 
general, Gg will involve a PV operator. 

To illustrate the PV formalism in these equations, we now specialize to mo- 
mentum space. Also, we assume that V is essentially a one-body potential, which 
can arise in several types of problems. For example, one may have a two-body 
problem that clearly separates into center-of-mass and relative coordinates, so that 
the potential is a function only of the relative coordinate. Another possibility arises 
when many-body effects can be absorbed into a self-consistent single-particle poten- 
tial. Moreover, we assume that the angular contributions can be removed from the 
formalism, with the wave function and operators decomposed into partial waves. 

With these specializations, Eq. (4.2) becomes®? 


PEO = plo) 4 Gio) yO, (4.4) 


where £ labels the angular momentum value, ko labels the energy 
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and giro) is the £-projected plane wave 


6(ko — k) 
ke’ 
with Ne(ko) a normalization factor. Note that the function on the right-hand side 


of Eq. (4.5) is essentially the nonangular part of the Fourier transform of a plane 
wave, which is a delta function,’ namely 


EY Yen) Yin(bo), (4.6) 


Lm 


(k | 6) = of" (k) = Nelko) (4.5) 


6(k — ko) = 


where Yem(k) = Yom(@x, Ox) is a spherical harmonic. Also, in Eq. (4.4) Vz = Ve(k, k’) 
is nonlocal in momentum space and Gh) is the Green’s function. 


P 


GY (k, k') = Tgm" —k’). (4.7) 

It is easily verified that®!-®° 
Vedi? = Re of, (4.8) 

where the R-matrix is given by 
Re = Ve + Ve GEP Ry, (4.9) 


or explicitly in matrix form 


P 


TET Va(k,k') Re(k',ko). (4-10) 


2 mr 
Re(k ko) = Vlk, ko) — Fe f dk 


In Appendix B we show how to solve Eq. (4.10) numerically.?9 
We can also define another Green’s function®)™: 


2 -1 
GH ER) = [FB -! +ie)| e- 


= git) _ sr in6(K3 — k?)5(k k’), (4.11) 
the last line of which follows from Eq. (2.5). We then have® 
x = A + oval, (ag 
vx = Traps?) (4.13) 
and 
Ti = Ve + Ve GE Ty, (4:14) 


in analogy to Eqs. (4.4), (4.8) and (4.9). 


The Mathematics of Principal Value Integrals and Applications 861 


From Eqs. (4.11) and (4.14), we find that the T matrix, Ty, is given by 


P 


r gj Vek, k') Te(k’, ko) 
ð 


Te(k, ko) = Velk, ko) — = a| 
0 


. m 
— ingzkoVe(k, kat) , (4.15) 
and, as before, we remove the PV, obtaining 
2m [~ dk 
Te(k, ko) = Ve(k, ko) — FS E) 


at k'? Volk, k') Te(k', ko) — k3 Volk, ko) Te(ko, ko)] 


E im T3 ko Vilk, ko) Tao. k)} f (4.16) 


In Appendix B, we also show how to determine the T-matrix in momentum space. 
For additional results involving the R and T matrices, see Refs. 28, 45 and 51. 


5. Applications to Nuclear Physics 


We now discuss some results from intermediate-energy nuclear physics, in which the 
loop integral in Fig. 1 is evaluated for nucleon-delta-pion interactions (see Refs. 12, 
13, 16, 52 and 60). The aim is to solve iteratively Dyson’s coupled equations for 
this system. There are two-loop integrals that are of great interest in the theory: 
(a) the nucleon-delta integral, which becomes a self-energy insertion into the pion 
Green’s function or propagator, and (b) the nucleon-pion integral, which becomes 
a self-energy insertion into the delta particle Green’s function. 

In our model, we will assume that we are dealing with infinite nuclear matter at 
zero temperature with all propagators assumed to be causal. Also we will assume 
that the nucleon propagator, Gy, is taken from a mean-field model, with 


O(en(p) — un) ih O(un — en(p)) 


(w—en(p) +e) (w—en(p) — te)’ (5.1) 


Gy (p, w ) = 
where p is the spatial momentum, with p = |p|, and w is the frequency (or energy); 
Ln is the nucleon chemical potential. In Eq. (5.1), on the left-hand side we indicate 
a general dependence on p, while the right-hand side in actuality involves p. This 
behavior will be true of all nuclear matter functions discussed in this section.1® In 
this theory, we will also assume that the pion and delta propagators can change 
with iteration, but the nucleon propagator will always be specified by Eq. (5.1). 

The mean-field energy in Eq. (5.1) is given by 


h2 


= —— +V (p) +m e, 5.2 


en(p) 
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where my is the nucleon mass, and the effective nucleon mass, m% , and mean-field 
potential, v, both depend on the nucleon matter density 


2 3 
p(pF) AT, 372 PF, (5.3) 
with pp the Fermi momentum. In Eq. (5.1), O(z) is the usual step function, while 
£ is infinitesimal, so that from Eq. (2.5) we obtain 


(w ~ ew(p) + ie)" F 6(w — en(p)). (5.4) 


Zye 
e—=0+t w — en (p) 


Reference 16 presents a method for calculating mẹ and vO, using a Skyrme-force 
model. 

On the first iteration, we assume that the imaginary part of the delta particle 
propagator is given by!?16 
1 Ta(w)O(w — us) 
2 [(w — ea(p))? + TA)’ 


where pa is the delta chemical potential and e, is the delta mean-field energy 


Im [Ga(p,w)] = - (5.5) 


2 
eal) = FG + VA) +mae, (5.6) 
with 
ma(p) _ miy(P) (6.7) 
ma MN 
and 
VL (p) = VP(p). (5.8) 


The restriction in Eq. (5.5) tow > wa is to prevent the propagation of holes. 
Reference 16 gives a prescription for calculating the delta width, [a (w). Then, we 
calculate the real part of the delta propagation using Eqs. (3.1)-(3.3), giving 


Re [G1 (p, w)] = 2 [ dus! MCap, w’) sign(w" = Ha) (5.9) 


(v =w) 


which can also be evaluated in closed form.!® Equation (5.9) rigorously guarantees 
that causality is satisfied for the delta particle. 
We can evaluate Eq. (3.36) in the limit of zero temperature,!?"!® to obtain 


Imiéva(pw))=-2 fd! fea 
x [O(w +w — pn) O(pa — w) + Ofun -w — w')O(w’ — pa) 
x [F(|2q + p|, 2w’ + w)|?Im[Gy(p + q,w + w’)|Im[Ga(q,w’)], 
(5.10) 
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where we have also made two generalizations: (a) including a three-dimensional 
spatial momentum integral that is “covariant” with the frequency integral and (b) 
inserting the form factor 


= ; A, B consts, 


F(p, w) = (w? + \p|?)? + B’ 


in order to assure convergence of the analogous nucleon—pion loop integral (to be 
discussed). We emphasize that in Eq. (5.10) the form factor is not really crucial for 
convergence. However, it is necessary in the nucleon—pion integral, for consistency, 
it can also be used in the nucleon-delta integral. 

We recall that the nucleon—delta loop integral is to be inserted as a self-energy 
correction into the pion propagator. The actual insertion involves several additional 
manipulations.!®5? First, we must define 


Im [UPA (p,)] = Im [Cva(p,w)] + Im [Cva(-p,w)] - (5.11) 
Then, from Eqs. (3.38)—(3.41) we find, in the zero temperature limit, that 


J=2 fa im [Upao] stent’) 


0 
Im [UPA (P ) (ew) 


(5.12) 


in which we have also assumed that 


HA ZHN, 


which appears to be satisfied for densities up to about twice normal nuclear matter 
density.! As before, Eq. (5.12) guarantees that causality is rigorously satisfied for 


the pion self-energy correction. Next, the full vertex correction is given by!®*? 
Una (p, w) = [1 + sUy; dip.w)] UŞA (P w), (5.13) 
where 
iz 8r? r 
O(he}miypr ^ 


and gi, = 0.45 is the Migdal parameter.!® In practice, it has been found that Uya 
and uv) differ negligibly.1® The largest difference is ~ 10 % for p < 1.0 fm™, with 
the discrepancy decreasing with increasing p. 

Finally, the explicit expression the pion propagator, including the self-energy 
insertion from the nucleon-delta loop integral, is 


D(p,w) = (1 T Dalp,w)Er(p, w) t Do(p,w) a) (5.14) 
where 


Do(p,w) = (W? — pe? — m2c* + ie)! (5.15) 


864 K. T. R. Davies et al. 


is the usual free-particle casual pion propagation (with ¢ infinitesimal), and 


gava 
Bpa) = -Eho (F) Uva(pwhs (PP=1287, 610) 


is a “polarization” operation.!ê 


After a series of additional manipulations, a pion-nucleon loop integral, PIN 
(p,w), can be defined in analogy with Eq. (5.10). However, in this integral the 
pion propagator is multiplied by an additional function involving Una (p, w). Also, 
as mentioned before, for this loop integral the form factor is crucial for obtaining 
convergence. Then, after the imaginary part of Ua(p,w) is computed, we obtain 
its real part from 


(0) -P [7 p m EaP o’) sigaw — uw) 


and = (p w) is inserted into an iterated delta propagator!®: 
[Ga (p u)? = w — ea(p) - D4) (P w). (5.18) 


Thus, >) (p,w) is the exact self-energy correction to the delta propagator. Then, 
in principle the entire process outlined above can be repeated, until we achieve some 
type of self-consistency. Here we just deal with the first iteration. The imaginary 
and real parts of Uy a (p, w) versus w are plotted in Figs. 3 and 4, respectively, for five 
different p values.! The curves Re(Una)’s have the expected kind of behavior for 
functions computed from the dispersion relation, Eq. (5.12). That is, for functions 
behaving as in Fig. 3, their dispersion relation “partners” have the types of shapes 
displayed in Fig. 4. 


IMAG [Unga (P.c0)} x 10° (MeV?) 


10.0 
w x 10°2 (MeV) 


Fig. 3. The Im [Una(p,w)] vs w for p = 0,2,4,6 and 8 fm? (the full, dotted, dot-dashed, 
short-dashed, and long-dashed curves, respectively). 
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REAL [Un, (p.co)] x 10°3 (MeV2) 


“0.0 5.0 10.0 15.0 20.0 
@ x 10? (MeV) 


Fig. 4. The Re [Una(p,w)] vs w, for various p values. The conventions for the curves are the 
same as in Fig. 3. 


IMAG [E'2)ip,ca)] x 10°3 (MeV) 
b ; 
> 


-0.5 


w x 10°? (MeV) 


Fig. 5. The Im [EP pw] vs w, for p = 1,3,5,7 and 9 fm! (the full, dotted, dot-dashed, 


short-dashed, and long-dashed curves, respectively). 


The imaginary and real parts of D0) (p,w) are plotted in Figs. 5 and 6 (see 
Ref. 16). Again, the dispersion relation partners displayed in Figs. 5 and 6 have the 
expected shapes. However, we note that the behavior with respect to p shown in 
Figs. 5 and 6 is very different from that shown in Figs. 3 and 4. 
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0.4 


0.2 


0.0 


REAL [£'9)(p,ca)] x 10°3 (MeV) 


-0.4 


-0.6 
0.0 


@ x 10? (MeV) 


Fig. 6. The Re [ELP w)| vs w, for several p values. The conventions for the curves are the 
same as in Fig. 5. 


In any case, here we are not trying to resolve the complicated questions in the 
theory of pion-nucleon-delta-interactions. Our main purpose is to show that the 
formalism of two-particle loop integrals, discussed in Sec. 3.3, has applications to 
real nuclear physics problems. In turn, this formalism can be derived, as shown 
previously, using the PB theorem, which shows once again that the theorem is 
applicable to many types of problems in physics. 


6. Applications to Transport Theory 


Transport equations describe the rate of change of macroscopic physical quanti- 
ties (e.g. number of particles, density, energy, electric charge) due to convection, 
leakage, absorption, etc. The fundamental transport equation was postulated by 
Boltzmann in 1872 to describe the temporal evolution of the one-particle distri- 
bution function of a rarefied gas with binary collisions.? Thenceforth the nonlinear 
integro-differential Boltzmann equation has become the paradigm of transport equa- 
tions and has been adapted, simplified, generalized, or extended to cover various 
situations. For instance, linear or linearized approximations are used for neutron 
transport and gases close to equilibrium. Grazing collisions that occur in electron 
scattering lead to small angle expansions and the ensuing Fokker—Planck type equa- 
tions. Applications to charged particles introduce nonlinearities in the convection 
term itself (as opposed to the collision term) and are treated in the Vlasov equa- 
tion formalism. Inclusion of several species of particles form systems of coupled 
Boltzmann equations, finite size effects are taken into account within the Enskog 
approximation, and higher order collisions that take place in denser gases can also 
be accommodated. The interested reader can find the basic principles of transport 
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equations and their developments in one of the excellent monographs in this field 
(see Refs. 8, 9 and 19, and references therein). 

It is beyond the scope of this section to review all the situations described by 
transport theory whereby the PB formula or its generalizations can be applied. 
Instead, we shall restrict ourselves to illustrate these applications to a couple of 
simple examples that also bear historical significance. Since the collisionless Vlasov 
equation and the monoenergetic neutron transport were the first two transport 
problems where the PB formula has been recognized, (sometimes) rederived, and 
applied, we shall focus on these cases. 


6.1. Collisionless electron plasma 


In the framework of transport theory, the electron plasma is a charged gas described 
by the one-particle distribution function fo(v) + f(r,v,t). Here fo(v) is the equi- 
librium distribution function, which is assumed to depend only on the magnitude, 
v, of the velocity v, while f(r, v,t) describes the perturbation, depending on posi- 
tion r, velocity v, and time t. The ansatz fo = fo(v) is equivalent to saying that 
equilibrium state is homogeneous and isotropic. Assuming that the plasma is very 
rarefied (i.e. it is in the collisionless or Vlasov regime) and linearizing around fo, 


the transport equation for the electron gas can be written as® 33; 
of of e ô fo Re 
So oe m w ea 
E — r— r’ d 1 U 1 d 1 
(r,t) =e jr-r®| r f(r V ,t) V, (6.2) 


where e denotes the electron charge and m its mass. Looking for solutions of the 
form: 


f(r,v,t) ~ eket drul), (6.3) 
we obtain 
4re? k- ð ! 1 
(w -k v)órw(y) = -2 Ih Sarul )dv’. (6.4) 


Equation (6.4) can be further simplified. We take k as the z-axis, denote —¥ = 
vz, = u, integrate over vz, vy, and use the identity 


f 7 ae f k do -$h zh (6.5) 
to find 
(Wol En S Broi, (6.6) 
where ~ 


Prw(U) = T dv, ia duy Pkw (V) (6.7) 
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and 


8r?e? u 


zfo(u) f (6.8) 


n(u) = = 


The integral f a y.,(u")du’ cannot be zero; indeed, assuming that this quantity 
vanishes, we find that ¢(u) = 6(w — ku) and thus 


a Pry (udu É 0, 


yielding a contradiction. Thus we may normalize the integral to one. With this 
choice, there are two classes of solutions® °° for Eq. (6.6): 


Class 1. 
A ee (6.9) 
A w— ku’ 
for all the discrete w that satisfy the equation: 
oO feu) 
Ay (w) =1 Je ae ba te =0. (6.10) 
The function A;,(w) is usually called the plasma dispersion function. 
Class 2. 
- _ [p mlu) _ 
Pkw (u) = le Tku + A(k,w)6(w — ku)| . (6.11) 


The arbitrariness in choosing À(k,w) can be utilized to get a complete set of func- 
tions and to satisfy the normalization constraints. The latter implies that Pk w(t) 
does not contain other singular terms — for instance derivatives of the 6-function. 
The normalization condition f% 4, „(u)du = 1 yields A(k,w) =1-P [S du. 

The generalized eigenfunctions (6.11) are complete and satisfy orthogonality 
relations in the form’-°; 


J Pkw lU)Pr w (U)ny (u)du = N(k,w)d(w—w’'). (6.12) 
To compute the normalization factor N(k,w), one applies the PB formula which 
yields’: 


d?(k,w) + 12172 (w/k) 
ne (w/k) 


We shall give more details of the application of the PB formula in Sec. 6.2. 


N(k,w) = (6.13) 


6.2. Monoenergetic neutron transport 


The neutron transport equation describes the advection, absorption, and production 
of neutrons in a host medium. The interest in this application grew dramatically 
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during and after WWII in relation with civilian and military applications (nuclear 
reactors and atomic bombs, respectively). 

The neutron gas is neutral and is assumed to collide only with the host medium. 
Thus, the corresponding transport equation is linear. Moreover, we shall assume, for 
simplicity, that the transport is stationary, the scattering is isotropic, the neutrons 
are monoenergetic, and the geometry is one-dimensional. Under these assumptions, 
the Boltzmann equation can be written®*4 


1 
ulije J Fæ d, (6.14) 


where f(x, u) is the stationary one-particle distribution function depending on posi- 
tion x and cosine y of the angle made by the velocity with the z-axis. In Eq. (6.14) 
the neutrons’ speed and the absorption cross-section are normalized to one. The 
constant c describes the multiplication properties of the medium (c < 1— subcriti- 
cal; c = 1— critical; c > 1— supercritical). We restrict ourselves to the case c < 1 
and look for solutions in the form: 


f(z, u) = e? p(n) , (6.15) 


which leads to 
ce f! 
w= we =F f bua (6.16) 


Since fi $ du cannot be zero (see the argument after Eq. (6.8)), we normalize it 
to one and obtain again two classes of solutions: 


Class 1. 
cy, 1 

bv, (MH) = 2 n = w (6.17) 

where the two discrete roots v; ¢ [—1, 1] are the solutions of the dispersion equation 
1 

_ cy 1 = 

Class 2. For v € [—1,1] we get 
cy 1 
(u) = Saran + A(v)6(v — u). (6.19) 


These generalized eigenfunctions satisfy the orthogonality relations 


1 
f h (Dhu (udu = N(Y) =v). (6.20) 


To obtain the normalization factor, N(v), for the continuum modes, we follow 
essentially the derivation in Ref. 8. The expansion of an arbitrary function f(u) in 
the continuum modes reads 


1 
flu) = [Alb war’ (6.21) 
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Multiplying by u, (u4) and integrating over u, we get 


Í wor (u) f(u)du = ie uo (p)du (f. A(v')64 (ud) (6.22) 


The left-hand side of Eq. (6.22) is defined to be the product between the expansion 
coefficient A(v), and the normalization factor N(v). The double integral on the 
right-hand side is evaluated by using the PB formula. Carrying this out, Eq. (6.22) 
becomes 


nae) = f nau |Z + ono - 0] 


«(fae eae Ti 


cv P 
= d mies a alae , 
[ame aS at 
i nea P 
+ k T amn- ; a 


z2 


Pf AAS = mar 
+ T D. A(v')A(v')6(v' — p)dv'. (6.23) 
i -1 


The last three terms in Eq. (6.23) are computed without difficulty, yielding: 


cv 1 P P 
Dwf. so zá )v' dv 
w f P 
-— —— A(v')A(v' )v'dv' + vA(v)A? (v). (6.24) 
2 ~1 y =p 


From the first term, by applying the PB formula (2.35) and identifying f(y, v’) with 
co pA(v'), we get 


1 1 1 
hae du- P “al dy! = cw Prien. 
—1 a 


P Êw 7 
—— ; .25 
S ay) + ia dv’ Pay carer ——pA(v')dp. (6.25) 


In (6.25) we make use of the identity 


(6.26) 


The Mathematics of Principal Value Integrals and Applications 871 


to get 


2 fl P 9 P P 

Me 1 1 1 _,' 

a faa [au pe ” =] 

= oof pce Zion — A(v))| v' A(v’)dv’ (6.27) 
~ 4 fy v-vite 


When collecting the results, the terms in (6.27) cancel the first two terms in (6.24). 
2 
Eventually, we obtain N(v)A(v) = vA(v)A?(v) + 7? A(v), which yields 


12 22 
Nw) =» [20) + - | . (6.28) 


A similar derivation was used to obtain (6.13). Other applications of the PB formula 
to transport theory can be found in Refs. 36 and 43. Related applications to the 
half-Hilbert and half-Hartley transforms have been recently discussed in Ref. 50. 


7. PV Integrals in Condensed Matter Physics 


Because the electron energy distribution in many solids is continuous in character, 
multidimensional PV integrals arise naturally in calculations relating to condensed 
matter physics. This is particularly the case for those of a perturbative nature, 
where vanishing energy denominators are the rule. Due to the need for high accuracy 
in distinguishing subtle features, the numerical treatment of these integrals has 
often strained the limits of computational technology. In spite of this, some results 
remained tainted with controversy until exact analytical evaluations were achieved. 
Even in discrete lattice problems what is interesting, in many cases, is the real part 
of a causal Green function, often best expressed as a PV integral over a polyhedral 
domain in reciprocal space. This section contains several examples of such principal 
value calculations important in solid state theory, and their analytical resolution. 
In each case, the interesting features are the means by which the singularity in the 
integrand is treated analytically and the intricate analytic structure revealed. 

Forty years ago, Lindhard?’ evaluated the density response function for the non- 
interacting homogeneous electron gas, Io (k, w) — one of the fundamental quantities 
in the theory of metals — represented by a Feynman diagram similar to the one in 
Fig. 1. 

A formal expression for the lowest order corrections, due to Coulomb interactions 
among the electrons, corresponds to the sum of three diagrams like the one in Fig. 1, 
two with Coulomb self-energy parts in one of the two propagators, and the third 
having a vertex correction connecting the two propagators. This was worked out 
for the first time by Dubois!® in 1959. In the static case (w = 0), in suitably scaled 
form, this is proportional to the six-dimensional PV integral 


872 K. T. R. Davies et al. 


Here, R is the region {k||k — $q| < 1}. The numerical study of I(q) began 
with Geldart and Taylor, who used the Pauli principle to isolate the singularity 
and were able to calculate the result in the range 0 < q < 2. Because of questions 
concerning gradient terms in density functional theory, the precise behavior at q = 0 
and q = 2 is required. In 1990 Engel and Vosko?’ derived the analytic expression 


1(q) = —*=#* (m ah tak [ols (mI) 


aese [alse (WB) 


The method used to obtain (7.2) was to decompose J(q), defined in (7.1), into 
four integrals, two with denominators (k-q)?(k+k’)*, which can be evaluated easily 
in terms of spherical coordinates, and two with denominators (k-q)(k’-q)(k+k’)?, 
which cannot. The first of these, A_(q), is even in q and reduces easily to a fourfold 
PV integral in cylindrical coordinates. Engel and Vosko then showed that it satisfies 
the differential equation 


DA_ = 


7.3 
tam (7.3) 
with boundary condition limg—oo q* A_(q) = const, where D is the differential op- 
erator 


Pau eee. (7.4) 


The unique solution to this problem is found directly. The attempt to carry 
this procedure out for the second integral, A;(g), failed due to the more complex 
q dependence of the integrand. To get around this, Engel and Vosko expanded 
A+(q) to order q~® for large q and noted that the resulting series is, to that order, 
consistent with the differential equation 


DA,(q) = 


1 1+ 
Fi rar (7.5) 


qi-@) 
for which the unique solution can again be found directly. Thus, while the deriva- 
tion of (7.2) relies on the unproven assumption (7.5), there is little doubt that it 
is indeed exact, since it reproduces all previously known information about I(q). 
Recently, (7.1) has been evaluated exactly?” by the reduction of (7.1) to a single 
integral which can be evaluated in terms of PV integrals of elementary functions 
and polylogarithms.*° The PV singularity in terms where it could not be eliminated 
was evaluated by means of the definition. Wherever possible, it was first reduced 
to a logarithmic (and therefore integrable) singularity through integration by parts 
producing integrals with one logarithm in the integrand, leading to the Euler dilog- 
arithm, and integrals with two logarithms, leading to the trilogarithm. The result, 


for 0 < k < 2, is 
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I(k) = ~(4—k?)In 


i zee) | ale- (235) 
ren (222) Lae =) i ( 25) 


— 41n(2) In =) a 5 (2— b) in? (Fz) 


24k 

+ alri (33) +r (27) +21n?(2) ~ 7/3} = saa! - k?) 

: {agit k?) In? G=) — 6(1 — In(2)) In? (=) 

= sle + k) ln? ==) — (2 — k) ln? Ga] 

m Te - b)[ In (7 z =) In (==) — in (=£) In (=) 

A aleran ==) + (2— k) n? E3) - 61n2(2)In(2=*) 
+ 


6 In(2)(2 — In(2)) In ( zË) 


+Žm (8), (48 eruan e- k)In(2 — k) — 4] 
— 6ln 24t) In (=) In (; z£) + 2] Ge) 
zii (e4) - a -m@) [tn (244) +1 (275) 
CE) 


+ 27? (1 — ln(2)) — 121n?(2) + 101n3(2) — xo. (7.6) 


By examining this formidable expression, one can appreciate the intricate analytic 
complexity lurking behind the PV symbol, which is frequently missed in numerical 
studies. By comparing (7.6) with (7.2), expressions for Engel and Vosko’s integrals 
were deduced and proven by symbolic differentiation, thereby proving the validity 
of the assumption on which (7.2) is based. 
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Let us now turn to the second class of problems which result in multidimensional 
PV integrals. Consider a regular lattice for which the nearest neighbors to a given 
lattice point, which we assume to be the origin, are 6;,...,6.. The corresponding 
lattice Green function is then the solution to the difference equation 


c 
zG(z;r) -J JjG(z;r + 6;) = bro, (7.7) 
j=l 
where z is a complex parameter and J; is a fixed “bond strength”. This function 
can be expressed formally as 


e 
G(z;r)= q Je Tk 
(7.8) 


w(k) = > vies, 
j=1 


where the integral is over the unit cell of the reciprocal lattice, having volume Q. 
These functions play a central role in crystal physics where one generally writes?” 
z = t — ic with t real and ¢ infinitesimal. The imaginary part 


1 
p(t) = = Im G(t — ic; 0) (7.9) 
plays the role of a density of states, while the real part 
ek r 
(tir) = È 5 fa ae 5 (7.10) 


acts as a generating function for transport paths for elementary excitations on the 
lattice. In particular, LOD J;,0) is directly related to the probability that a 
lattice excitation eventually returns to its starting point. The high precision evalu- 
ation of these PV integrals for lattices in three and higher dimensions is extremely 
difficult and is the subject of hundreds of studies since special cases for the three 
cubic lattices were first worked out exactly in a classic paper by G. N. Watson*® 
1939. Unfortunately, no comprehensive review of these studies is available, but the 
earlier work is surveyed and documented in Ref. 32. Here we shall look at exact 
expressions for the “simple cubic” lattices (J; = 1) in various dimensions. 

In one and two dimensions the PV can be found by analytic continuation. We 
have, in the one-dimensional case, 


1 kig 
t = — 
stn) == Re | 


Now, for t > 1 the integral is [t— Vt? — 1]"/Vt? — 1, so by analytic continuation 
to0<t <1, we find 


cos(nx) 


t — ie — cos(x) (oan 


_Im [t+ivi-#]" Un) 
vize — Vize 


g(t; n) = (7.12) 
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where U,,(t) is a Tchebyshev polynomial of the second kind. Note that the corre- 
sponding Watson integral (n = 0,t = 1) does not exist. 
In two dimensions, the Green function is 


G(t, (m,n) =a f° in __cos(mky) cos(nka) (7.13) 


t — ie — cos(k,) — cos( k2) 


and for t > 2 can be written as the convergent Bessel integral 


is e* J,,(2)Jn(x)dx = Gna = 2) 


1 1 + + 
mint miat 1+ md 2 1+ mi n 


x 4F3 (7.14) 


m+i n+1 m+n+4+1 + 


Here „F, denotes the generalized hypergeometric function.?! The analytic continu- 
ation to the region 0 < t < 2 can be found®® and oe 


(-1)” l+m Bin l-m-n -n itm m+n 
G(t;(m,n)) = F; a “a 
2 i l 1 t 
2 2 4 
2 (—1)™+1(m? es n?)t mete moe +2 amon ae 
4 3 1 t 
2 2 4 


(7.15) 


if m and n have the same parity, for example. Again, the corresponding Watson 
integral (t = 2,m = n = 0) does not exist. In three dimensions the matter is 
drastically more complicated, and only a few special cases have been worked out 
analytically. Let us look at the extended Watson integral 


Pf [E dk 
Mes af f f 1 — å(cos(k1) + cos(k2) + cos(k3))’ (7:16) 


for which the integrand is singular for —1 < z < 1. This was evaluated by Watson*® 
for z = +1 in 1939, but in spite of numerous attempts, was not evaluated for 
general z until 1973 in an important paper®? by Joyce. The following is a somewhat 
simplified version of his calculation. 

Formally, the pole in a principal value integral can be eliminated, at the cost of 
an additional integration, by the identity 


1 o0 
= f e °7e*ds, (7.17) 
v—y 0 


where the convergence, if x < y, depends on the nature of the original integrand. 
By applying this to (7.13) and invoking the integral representation for the modified 
Bessel function, we obtain 


W(z)= L e™° IÈ (sz/3)ds. (7.18) 
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Next, we use the identity?® 


S (1/2) 
B(s) =Y E 4 y[-k, —k, 1/2; 1,1; 4]a* (7.19) 
2~ (2k)Iz! 
to write W(z) = >, Anz?”. The coefficients A,, being hypergeometric, must 
obey a simple recursion relation,?® which in this case is 


36(n + 2)? Anya — 2(2n + 3)(10n? + 30n + 23)An+1 
+ (n+1)(4n? + 8n + 3)A, = 0. (7.20) 


Identifying (7.20) as the recurrence relation arising in the series solution of a differ- 
ential equation, we find that W(z) must be a solution to the third-order ordinary 
differential equation (z = z?) 


W” (2) + 3f(2)W" (x) + [2(f(x))? + f'(z) + 49(x)|W"(z) 
+ [4f(x)g(x) + 2g'(2)]|W (z) = 0, 


1 1 1 
f2)= 2+ 5G nh tao ey) 
3(x — 4) 


(2) = Tee — 1) =)" 


It is known from the theory of generalized hypergeometric equations that the so- 
lution to (7.21) regular at the origin is proportional to [y(x)]*, where y(x) is the 
regular solution to the second-order equation y” + f(x)y’ + g(xz)y = 0. Joyce was 
able to identify this with Heun’s equation, showing that W/(z) is essentially the 
square of a Heun function. Finally, by invoking various transformation properties 
of Heun functions,** he obtained 


(1 — ĉr) 4 


W(z) = al E 1) aK (ky) K(k) ’ 
1 
kį = 5+ isai — z2 — 702 —22)/1— 22, 
where 
1 1, 1 1 
= gta Sel (ie 32) o CETE 
m= 5 tee? -zyl - 2201 - 52%), 
Tı 
=; 7.21 
T2 zı -1 ( ) 


Once again, one sees the striking analytic complexity unfurling from the exact 
resolution of a simple looking multi-dimensional PV integral. 
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8. Summary 


In this paper, we have categorized and discussed various types of PV integrals 
that arise in applications. Both infinite-limit and finite-limit integrals have been 
examined, including comparisons for the former with complex variable theory and 
analyses for the latter using quadrature relations. For each of these categories we 
have considered the following cases: 

(a) one-dimensional integrals containing a simple pole. 

(b) one-dimensional integrals containing a single multiple pole. 

(c) double integrals containing two simple poles. 

(d) double integrals containing two multiple poles. 

(e) multiple integrals containing products of simple poles. 

Also, in many of these cases we have discussed or referred to various numerical 
methods for efficiently evaluating the integrals that occur. 

In these developments we have tried to clarify certain misconceptions regarding 
the PV formalism. First, we have shown that the PV can be considered as the 
convergent part of a divergent integral. Second, we have shown how to generalize the 
usual definitions applicable to a single pole in a one-dimensional integral to higher 
order poles and multiple integrals. The role and applications of the famous PB 
theorem have been discussed. Third, we have shown that the value of a finite-limit 
PV integral can be expressed as the difference between the quadratures evaluated 
at the end points of the range of integration. In this regard, the PV can be viewed 
as the value of the integral obtained by ignoring the singularity occurring in the 
integrand. 

Many important applications of the PV formalism have been reviewed. Various 
types of dispersion relations have been analyzed, including applications of the PB 
theorem to two- and three-particle loop integrals. The R- and T-matrix formalism 
in momentum space was briefly discussed. Then we have reviewed representative 
applications to nuclear physics, transport theory, and condensed matter physics. 

Finally, various important extensions of the work should be considered in the 
future. One interesting idea is to attempt to evaluate the convergent part of an 
integral containing singularities other than poles, e.g. branch cuts. Another impor- 
tant application of the work might be to connect the PV formalism with the Wick 
rotation?! in quantum field theory, which has been used for bound state equations 
in which Minkowski variables are analytically continued in order to generate a set 
of Euclidean coordinates. 
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Appendix A. Summary of Methods for the Numerical Evaluation of PV 
Integrals 


A number of methods have been developed for evaluating numerically a variety 
of PV integrals (see Refs. 10, 13, 15, 17, 29, 54, 55 and 60). Also Refs. 17, 54 
and 55 contain additional references to numerical methods. In this Appendix we 
summarize various numerical techniques used to evaluate PV integrals. 

(a) Single Integrals of the form: 


h(x) =P a aE (A.1) 


(2'2 — x?) a x?) 


References 10, 29 and 60 treat such integrals. In Ref. 10, a complex-plane 
method is used to evaluate integrals like (A.1) when f(x) is a highly oscillatory 
function. 

(b) Single Integrals of the form: 


Io(x) =P 3 f Ua. (A.2) 


oo Vr 


In Ref. 13, methods are presented for evaluating Eq. (A.2). Also, special algo- 
rithms are developed for the following cases: 


(i) when f(x) contains a single step function, namely 
f(x) = O(z — u)g(z) (A.3) 


and 


(ii) when f(z) contains two symmetric step functions: 
f(x) = [O(x — u) + O(—2 — p)} 9(z). (A.4) 


Applications involving zero-temperature causal Green’s functions for fermions and 
bosons furnish examples of Eqs. (A.3) and (A.4).13 (See Sec. 3.1, in particular 
Eq. (3.3).) 


(c) Single Integrals of the form: 


koer a (A.5) 


s'r’ 


where b and c are finite. Algorithms for evaluating Eq. (A.5) are given in Refs. 13 
and 15. 


(d) Single Integrals of the form: 


_ c f(2')da' 
Is(x) =P gaa (A.6) 
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where b and c are finite and n > 1. A method for treating this case was presented 
in Ref. 15. Because of the formalism developed in Sec. 2.1.3 for higher-order poles, 
there are boundary terms that involve derivatives of f(x). Thus, in order for this 
numerical method to work, either these derivatives must be obtained analytically 
or the function f(x) must be amenable to accurate numerical differentiation. 


(e) Double integrals of the form: 


P P 
Iu) = | ir f av fla) (A.7) 


where the integration limits can be either finite or infinite. The numerical evaluation 
of such integrals is still being perfected. Some methods for calculating integrals like 
Eq. (A.7) were presented in Refs. 15 and 60. Also, precision calculations of loop 
integrals (which involve two PV integrals, as discussed in Sec. 3.3) were discussed 
in Refs. 13, 14, 16 and 60. 

One of the challenges of the calculation of Eq. (A.7) is to compare it with the 
numerical evaluation of 


E N E E 
hw) = fav |g” te. (A.8) 


The difference between Eqs. (A.7) and (A.8) is the PB “interchange term” given 
in Eqs. (2.35) and (2.36). This interchange term vanishes if at least one of the 
PVs in Eqs. (A.7) or (A.8) does not occur (e.g., if the point x = u is outside the 
range of integration). This property has been nicely demonstrated for an analytic 
example,!® but it is instructive to exhibit it for cases in which the integrals must 
be evaluated numerically. 

All numerical methods for evaluating PV integrals make use of some type of 
subtraction trick in order to eliminate the formal PV of the integral. In this way, 
we rely on the property (discussed in Sec. 2.1.4) that the PV is the convergent part 
of the integral. We illustrate by evaluating Eq. (A.2) or (A.5), defining’? 


° f(x) dx 
(x —a)’ 


where a is real, b < c, and b and c are both finite or both infinite. Equation (A.7) 
can then be rewritten as 


I(a)=P (A.9) 


[f -i f(a) 
=P [& =a “t f(a) P [5 (A.10) 
Now the PV symbol can be removed in front of the first RIR so that 
I(a) = f fore dx + (aye f ay (A.11) 


If c > +00 and b > —oo, we have from Eq. (2.62) 


pf BN Li, (A.12) 
t-a 
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and if both c and b are finite 
© dr c—a 
= ln |—— A.1 
Pf z-a alga! (ALa) 


a result which follows from Eq. (2.21a). Thus, the second integral in Eq. (A.9) can 
always be evaluated analytically. 

In general, though, the first integral in Eq. (A.11) will have to be evaluated 
numerically. See Refs. 13 and 29 for a discussion of the types of Gaussian meshes 
recommended for many integrals of physical interest. (See also Eq. (B.5).) After 
discretizing, the first integral in Eq. (A.11) can be written as 


Oe HeD - FO) (A.14) 


(xj — a) 


where N, is the number of Gaussian points and w, is the weight appropriate to the 
jth mesh point.}3-29 

However, Eq. (A.14) cannot be evaluated if a is chosen as one of the mesh points. 
This is a major disadvantage if we wish J and f to lie on the same mesh.!*® This 
difficulty can be overcome if f(x) can be differentiated and evaluated at the mesh 
points (e.g. by a numerical procedure). For such cases, we find that 


Ng 


T(z) = D Fe- fe), + | Wi, (A.15) 


j=1, jži (z; — zi) 


=T; 


a result obtained using L’Hospital rule.!° Unfortunately, it is not always possible to 
differentiate f(x), particularly if f(z) is a function that is not known analytically 
and that can only be determined by some prior numerical procedure. 


Appendix B. The Haftel-Tabakin Method for Numerically Evaluating 
the R and T Matrices in Momentum Space 


Due to Eqs. (4.5) and (4.8), we need only concern ourselves with the matrix values, 
Re(k, ko). (That is the right-most matrix index is always ko.) Following Haftel and 
Tabakin,?? we remove the PV in the second line of Eq. (4.10) by a subtraction trick, 
using Eq. (2.64), obtaining 


Re(k, ko) = Ve(k, ko) 


foo) 1 
Z f gg [EVel KRR, ko) ~ RBVeL I, ko) Reh, ko) 
h? Jo (k'?-— kô) (B.1) 
The R-matrix may be evaluated by a matrix inversion!®?%60 if Eq. (B.1) is 
discretized. In particular, we have 


No+1 
Vilki, kn 41) = D FIP (ki, k;)Relkj,kn,+1), (B.2) 


j=1 
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where N, is the order of Gaussian integration, kv, +1 = ko, and Fy is an Ng+1 X Ngi1 
matrix defined by 


FS?) (kj, ky) = ôi +O Velki, ky) (B.3) 
and 
E wkk — k2), for j < Ny. 
(R) _ 
Q; ai _ 2m (B.4) 


Ng 
KAY [wi/(k? — kj)] for j = N, +1. 


h i=l 
The Gaussian weight for point j is given by w;, and the quadrature used is extracted 


from the transformation!?:?9 
k; = tan F (uj +1), 1<j<N,, (B.5) 


where u; is the usual Gaussian—Legendre mesh point that lies in the range [—1, +1]. 
The mesh defined by Eq. (B.5) has been widely used to evaluate many types of 
integrals, particularly PV integrals. Overall, it is more accurate and efficient than 
other quadrature schemes used for evaluating infinite-limit integrals. 

By inverting Eq. (B.2), we obtain 


Ngt1 
Ra(kisky,+1) = >> FE" (i, kj) Velki, ky, 41)) (B.6) 


j-l 
which gives the R-matrix, including the on-shell case 
Re (kn, +1, Nn, +1) = R(ko, ko). 


Clearly, in this formulation, kg cannot be one of the Gaussian mesh points, which 
may be a disadvantage in some applications. See Ref. 60 for further discussions of 
this problem and an alternative, approximate formulation of the method. Also, see 
the discussion regarding Eq. (A.14) in Appendix A. 

Finally, we remark that, by comparing Eqs. (4.10) and (4.15), we can solve for 
the T-matrix in a similar fashion. We redefine the 2 in Eq. (B.4) as 


2m ; 
qa vkj (ky — kô) for j < Ng, 
a”) = mm e mri (a 
=TM, [w:/(k? — k2)] + =a ko for j = N, +1, 
i=1 
which gives, in Eq. (B.3), a matrix FP. Then, we obtain 
Nyt ee 
Te(kis ksi) = XO (FE) (his key ely, v2) (B8) 


j=1 
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Note that, if the potential V; is real, the R-matrix is also real, but the T-matrix 
is complex. This implies that in Sec. 4 the pro) functions can be chosen real, while 
the x functions are always complex. Also, the beauty of this method is that one 
obtains Re(k, ko) and T;(k, ko) for a fixed mesh of k values and, as a bonus, the 
“on-shell” value for k = ko. 


Appendix C. A Useful PV Integral Transformation 
Theorem. Leta; >0, bj (j =1,...,N—1) be any real numbers and 
N-1 


o(2) =2- E 


z- bj 
j=1 i 


Then for any function F for which the integrals exist, and y £ b; for any j, 


_p [A EBn p [? EO y 
rap f. z-y sii rea 


Proof. The equation u = ¢(z) can be written in the form 


N-1 
G(x)/ T] (x - 8) =0, 


j=l 


where 
N-1 


N 
G(x) =u J] (x -8;) — f(x) = [lz - ztu) 
j=l j= 


and f is a polynomial. Here G(x) has the distinct roots z;(u) which are differen- 
tiable with respect to u for —oo < u < oo. Now, for simplicity assuming that y 
coincides with no 6; and omitting the PV symbol on the appropriate integral, 


t= fit [Pent fe Ae a. 


N-1 


Changing the variable of integration to x = 2;(u) over the range b;_1 < x < bj, we 


find a 
i x; (u) 
=P f F(u 1 du. 
"FO 


But 


zi; d a 1 
Doa Sgu n= u — (y) 


This completes the proof. 
As a simple corollary, we find that for 


N-1 


v@)= 


j=l 
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als a ls ars 


As long as w(x), (x) are convergent, N can be infinite. 
As an example, take f(x) = s(x? +a”*)~!. Then for ọ(x) = x — 1/z, we find 


from 
99 x dx Ta 
P aLaaa aLa? 
-œ T +atr-y y +a 


an _ av? -1) dr ma? 
wo t+ (a? — 2)? +1- y ~ yf + (a2 — 2)y2 +1 


which remains valid for y = 0. 


that 


For 
1&1 1 
Ra 2 k OE kr) t * e+ (ery 
we get 
p / Se cot(i/z) dem __ Ta __ 
-œ [z —cot(1/z)|? +a?z-y a [y—cot(1/y)]}? +a?" 
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